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ABSTAACT  (Continued) 


pyroclastic  blast  flow  became  buoyant  by  mixing  and  heating  surrounding 
air.  Observations  indicate  that  much  of  the  flow,  covering  an  area  of  600 
km  ,  became  buoyant  within  a  time-span  of  1.5  min  and  abruptly  ascended  to 
form  the  giant  cloud.  Calculations  are  presented  for  the  amount  of  air  that 
needs  to  be  entrained  into  the  flow  to  make  it  buoyant.  Assuming  an  initial 
temperature  of  450°C  and  a  magmatic  origin  for  the  explosion,  these 
calculations  indicate  that  the  flow  became  buoyant  when  its  temperature  was 
approximately  150°C  and  the  flow  consisted  of  a  mixture  of  3.25  x  10**  kg  of 
pyroclasts  and  5.0  x  10*  kg  of  air.  If  sedimentation  is  considered,  these 
figures  reduce  to  l.l  x  10**  kg  of  pyroclasts  and  1.0  x  10**  kg  of  air. 
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I.  INTRODUCTION 


The  initial  paroxysmal  explosive  eruption  of  Mount  St.  Helens  on  May  18, 
1980  involved  a  complicated  sequence  of  landslides,  vertical  explosions  and  a 
devastating  horizontal  blast,  or  pyroclastic  flow.  All  these  events  took 
place  within  a  few  minutes  and  were  immediately  followed  by  the  rapid  ascent 
of  an  enormous  eruption  cloud  that  reached  a  height  of  nearly  30  km  within  20 
rain.  Data  from  military  satellites  allow  the  detailed  behavior  of  this 
eruption  column  to  be  studied.  This  report  presents  these  data,  considers  how 
the  giant  cloud  relates  to  the  initial  events  of  the  eruption,  and  interprets 
the  data  in  terms  of  theoretical  models  of  atmospheric  convection.  The  most 
significant  new  observation  is  that  the  plume  spread  laterally  at  spectacular 
rates,  forming  a  giant  umbrella  cloud.  The  umbrella  cloud  spread  15  km 
upwind.  The  implications  of  this  behavior  for  understanding  large-magnitude 
explosive  eruptions  are  discussed. 


II.  OBSERVATIONS 


The  magnitude  5+  earthquake  that  occurred  at  Mount  St.  Helens  at  15:32.2 
UT  on  18  May  1980  was  immediately  followed  by  landslide  failure  of  the  north 
flank  and  summit  region  of  the  volcano.  Military  and  weather  satellite  data, 
together  with  ground  photography,  have  enabled  detailed  documentation  of  the 
events  in  the  first  few  minutes  of  the  eruption  as  described  by  Rice  (1981) 
and  Moore  and  Rice  (1984).  Explosions  were  first  observed  at  15:32.7  from  the 
summit  region  and  from  the  moving  landslide  that  fed  a  directed  blast  flow  to 
the  north.  There  is  currently  a  controversy  over  whether  the  directed  blast 
should  be  described  as  a  surge  or  pyroclastic  flow  (Hobblitt  and  Miller,  1984; 
Walker  and  McBroome,  1983,  1984;  Waitt,  1984).  For  simplicity  we  will  refer 
to  the  phenomenon  as  a  blast  flow.  The  blast  flow  moved  at  more  than  90  m/s, 
overtaking  the  avalanche  at  15:33.8  UT.  The  blast  flow  reached  its  furthest 
point  between  15:36  and  15:37,  covering  an  area  of  over  600  km  .  At  about 
this  time  a  huge  convecting  column  was  observed  to  rise.  In  about  ten  minutes 
the  cloud  grew  to  form  a  giant  mushroom- shaped  cloud  (Figs.  1  and  3)  with  a 
stalk  diameter  of  over  20  km,  a  height  of  25  km,  and  a  mushroom  cap  of  70  km 
diameter.  There  is  good  agreement  between  the  width  of  the  cloud  (interpreted 
from  satellite  data)  and  the  position  of  the  outer  edges  (deduced  from  the 
ground  photographs)  (see  Fig.  2). 

At  least  four  distinct  vertically  rising  clouds  can  be  distinguished  in 
the  first  15  rain  of  the  eruption;  these  are  identified  in  Fig.  4,  which 
records  data  on  their  variation  of  height  with  time.  Cloud  I  represents  a 
coalescence  of  material  formed  by  several  explosions  that  occurred  as  the 
landslide  blocks  moved  away  from  the  volcanic  cone.  Photographs  from  Mount 
Adams,  53  km  east  of  Mount  St.  Helens  (Fig.  10.3  in  Moore  and  Rice,  1984), 
show  that  the  major  effects  in  this  period  involved  the  lateral  spreading  of 
the  blast  flow,  but  that  vertical  motion  also  occurred,  with  the  apex  of  an 
inverted  cone-shaped  cloud  rising  above  an  area  about  4  km  north  of  the  summit 
region.  The  rise  of  the  apex  of  this  expanding  region  is  recorded  as  Cloud  I 
and  had  a  mean  vertical  velocity  of  25  m/s. 


Photographic  mosaics  of  giant  cloud  taken  from  (a)  U.S.  Forest 
Service  (1981),  and  (b)  Rocky  Kolberg.  Both  views  are  from 
north-west  of  Mount  St.  Helens.  Kolberg  photograph  taken  at 
about  15:52  GMT.  U.S.  Forest  Service  photograph  taken  at  about 
15:54  GMT.  Both  clouds  subtend  a  horizontal  angle  of  about 
90°,  as  shown  in  Fig.  2.  Photograph  in  Fig.  lb  assembled  by 
James  F.  Kolberg,  President  of  the  Toledo  Poster  Company. 


Camera  locations  for  the  U.S.  Forest  Service  and  Kolberg  photo¬ 
graphic  mosaics =  Rays  to  the  outer  edge  of  the  mushroom  cap 
are  shown  for  each  photograph.  The  diameters  of  the  cloud  from 
interpretation  of  satellite  images  at  15:50  and  15:55  GMT  are 
shown ,  which  bracket  the  time  when  the  photographs  were  taken 
(see  Fig.  5). 
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Fig.  3.  GOES-West  satellite  image  of  giant  cloud  at  16:15  GMT.  White 
dots  give  outline  of  Washington  State.  White  fringe  to  cloud 
is  interpreted  to  be  a  consequence  of  condensation. 
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TIME  AFTER  15.00  (  minutes  ) 


Altitude  of  clouds  above  sea  level  versus  GMT  for  initial  stage 
of  May  18  eruption.  Four  different  clouds  can  be  distinguished 
(see  text).  The  horizontal  bar  represents  the  duration  of  the 
blast  flow  at  ground  level  (2  km  above  sea  level). 


Photographic  evidence  (Moore  and  Rice,  1984)  shows  that  a  second  cloud 
(Cloud  II)  began  to  rise  rapidly  from  an  area  12  to  13  km  north  of  the  summit 
above  Spirit  Lake  and  somewhat  north  of  the  Toutle  River.  Moore  and  Rice 
(1984)  interpreted  this  cloud  as  the  consequence  of  the  second  landslide, 
containing  parts  of  the  dacite  cryptodome,  moving  into  the  Toutle  River  and 
Spirit  Lake  area.  They  envisage  that  the  Cloud  II  explosions  were  produced 
when  the  landslide  block  slammed  into  the  north  wall  of  the  Toutle  River 
valley  or  shattered  and  interacted  with  water,  triggering  a  large 
hydrovolcanic  explosion.  The  only  two  data  points  on  this  cloud  record  a 
vertical  velocity  of  70  m/s  (Fig.  4). 

New  information  from  U.S.  military  satellites  has  yielded  unique  data  on 
the  ascent  and  expansion  of  the  main  mushroom  cloud.  Figure  4  records  the 
variation  of  height  with  time  for  the  column.  Two  episodes  of  growth  could  be 
distinguished  and  are  designated  Clouds  III  and  IV.  Figure  5  shows  the  out¬ 
line  of  the  cloud  as  observed  from  Earth-orbiting  satellites.  The  diagram 
shows  cloud  outlines  at  5-min  Intervals  and  gives  the  first  data  on  the  radial 
expansion  of  a  major  eruption  cloud.  An  example  of  the  satellite  view  that 
can  be  obtained  is  shown  in  Fig.  3,  which  shows  the  giant  cloud  as  observed 
from  a  GOES  geostationary  weather  satellite  at  16:15  GMT. 

The  initial  vertical  velocity  of  Cloud  III  Is  approximately  110  m/s  and 
is  clearly  much  more  energetic  than  in  the  earlier  explosions.  Inspection  of 
photographs  and  Fig.  4  Indicates  that  the  center  of  growth  of  the  cloud  was 
displaced  to  the  north  of  Mount  St.  Helens,  in  the  same  region  as  the  explo¬ 
sions  that  generated  Cloud  II.  However,  the  data  in  Fig.  4  show  that  Cloud 
III  cannot  be  related  to  the  explosion  that  generated  Cloud  II,  which  evident¬ 
ly  decelerated  and  could  not  have  risen  much  above  10  km.  There  is  a  three- 
minute  delay,  apparent  in  the  diagram,  between  the  early  explosions  and  the 
development  of  the  major  eruption  column  (Fig.  4).  The  blast  cloud  had 
reached  its  maximum  extent  between  15:36  and  15:37.5,  at  precisely  the  time 
the  major  cloud  began  ft3  ascent.  Descriptions  of  the  blast  (Moore  and 
Sisson,  1981;  Waitt,  1981)  and  abrupt  termination  of  the  zone  of  devastation 
(Kieffer,  1981)  indicate  that  the  blast  cloud  terminated  lateral  motion  as  it 
lifted  off  the  ground  and  began  to  rise  vertically.  Data  on  flow  velocities 
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Fig.  5.  Outline  of  the  cloud  at  5-min  intervals  as  observed  from 

geostationary  satellites.  Cross  is  estimate  of  initial  center 
of  cloud  ascent.  The  stippled  area  shows  extent  of  blast 
deposit.  In  the  upper  right  hand  corner  the  wind  direction  is 
shown  as  measured  at  Spokane  on  May  18  at  elevations  of  5,  8, 
and  20  km. 


(Moore  and  Rice,  1984)  show  that  the  flow  still  had  a  velocity  of  about  100 
m/s  as  It  approached  the  distal  edge,  with  little  indication  of  any 
significant  deceleration.  The  outer  limit  was  thus  determined  not  only  by  the 
flow  slowing  down  but  by  the  flow  becoming  buoyant.  We  conclude  that  the  flow 
ascended  to  form  a  cloud  when  its  density  had  become  lower  than  the  overlying 
atmosphere.  We  therefore  interpret  Cloud  III  to  have  formed  by  the  buoyant 
ascent  of  the  whole  blast  cloud  over  the  entire  600-km  area.  The  main  cloud 
was  not  in  the  strict  sense  generated  by  explosions,  but  by  decompression  of 
the  flow,  by  mixing  of  air  with  the  blast  flow,  and  by  sedimentation  of 
pyroclasts. 

The  origin  of  Cloud  IV  is  more  problematic,  since  it  must  have  originated 
several  minutes  after  termination  of  lateral  blast  movement.  At  about  15:43 
satellite  sensors  detected  a  pulse  of  new  hot  material  rising  vertically  above 
the  vent  (Rice,  1981).  Moore  and  Rice  (1984)  interpreted  this  as  the  onset  of 
the  plinian  eruption,  discharging  predominantly  high-temperature  juvenile 
magma,  which  then  continued  for  a  further  nine  hours.  Cloud  IV  extended  the 
initial  column  to  nearly  30  km  altitude,  but  subsequent  plinian  activity  only 
produced  a  column  between  14  and  18  km  high.  The  first  stage  of  plinian 
activity,  if  manifested  by  Cloud  IV,  may  have  involved  a  much  greater  magma 

C 

discharge  rate  than  the  subsequent  plinian  activity  (approximately  10J  m  /s, 

*1 

as  compared  with  6000  nr/s  according  to  Sarna-Woj cicki  et  al.,  1981). 
Alternatively,  the  initial  plinian  activity  involved  discharge  into  the  giant 
cloud  itself,  which  may  have  enabled  the  plinian  column  to  rise  to  an 
unusually  great  height.  Shortly  after  16:00  the  column  height  decreased 
dramatically  to  about  14  km. 

The  satellite  views  of  the  cloud  (Fig.  5)  show  that  the  mushroom  cloud 
developed  very  rapidly,  with  an  initial  radial  velocity  between  15:40  and 
15:45  of  55  m/s.  The  cloud  had  expanded  15  km  upwind  of  the  western  edge  of 
the  zone  of  devastation  by  15:50,  at  which  time  the  cloud  reached  a  stagnation 
point  upwind,  but  continued  to  grow  in  downwind  and  crosswind  directions.  The 
relationships  at  the  stagnation  point  are  remarkably  simple.  The  upwind 
stagnation  point  represents  the  distance  at  which  there  was  a  balance  between 
the  radial  expansion  velocity  and  the  local  wind  velocity.  N0AA  satellite 


data  (Sarna-Woj cicki  et  al.,  1981)  show  that  the  feather-edge  of  the  cloud 
extended  between  12  and  16  km  altitude.  The  wind  velocities  at  Spokane 
(Washington)  recorded  by  the  U. S.  National  meteorological  service  at  these 
altitudes  on  May  18  varied  between  12  and  33  m/s,  averaging  22  m/s.  Between 
15:50  and  16:10  the  horizontal  radial  expansion  velocity  of  the  cloud  was 
fairly  steady,  between  25  and  20  m/s.  For  the  next  hour  the  western  edge 
moved  eastwards  at  about  5  m/s,  as  the  eruption  weakened. 

One  effect  of  the  wind  on  the  cloud  is  to  gradually  push  the  center  of 
symmetry  of  the  cloud  downwind.  Figure  6  shows  the  distance,  at  each 
successive  time,  of  the  center  from  the  cross  marked  on  Fig.  5.  The  cross  is 
assumed  to  represent  the  center  of  symmetry  at  15:37.5,  when  cloud  ascent 
initiated.  After  15:45  the  center  of  symmetry  is  moved  downwind  at  a  fairly 
constant  velocity  of  14.8  m/s.  This  velocity  is  approximately  half  of  the 
observed  stratospheric  wind  speed.  The  easternmost  (downwind)  edge  of  cloud 
reached  a  steady  velocity  of  30  m/s  after  15:50  GMT;  this  cloud  edge  is 
thought  to  be  purely  wind-driven,  corresponding  to  the  maximum  stratospheric 
velocity  observed  at  10  to  12  km  on  May  18. 

Wind  direction  information  measured  at  Spokane  at  17:20  GMT  is  shown  in 
Fig.  5  at  elevations  of  5,  8,  and  20  km.  Initially  the  outer  edge  of  the 
cloud  (at  15:45,  for  example)  is  elongated  ENE  under  the  influence  of  the  low- 
level  winds,  but  at  later  times  is  elongated  E  to  ESE  as  the  winds  at  higher 
altitude  influence  the  cloud's  shape. 

A  prominent  feature  of  the  mushroom  cloud  can  be  noted  in  cross-sections 
of  the  cloud  seen  via  radar  observations  (Sarna-Woj cicki  et  al.,  1982). 
Although  the  cloud  immediately  above  the  source  region  reached  25  to  30  km, 
the  altitude  of  the  upper  surface  decreases  radially  away  to  about  20  km. 

This  behavior  is  typical  of  a  plume  intruded  into  a  stratified  environment 
(Turner,  1979).  The  momentum  of  the  cloud  takes  dense  material  to  a  height 
where  it  is  now  surrounded  by  lower-density  air.  The  material  then  flows 
outwards  and  downwards  to  Its  own  density  level.  Higher-altitude  winds  were 
also  relatively  light. 


Figure  1  shows  that  the  ascending  column  beneath  the  umbrella  region 
converges  inwards.  Eruption  columns  originating  from  a  localized  vent  expand 
with  height  as  the  ascending  plume  entrains  air.  Such  columns  can  be  approxi¬ 
mately  treated  as  if  they  emerge  from  a  point  source.  The  giant  cloud,  how¬ 
ever,  originated  from  a  large  area.  Turner  (1979)  has  discussed  how  inward 
flow  is  typical  of  buoyancy  generated  from  a  heated  surface,  because  a  thermal 
boundary  layer  is  created.  The  fluid  streamlines  in  this  boundary  flow  con¬ 
verge  into  an  ascending  region  that  is  narrow  compared  to  the  width  of  the 
source  area.  The  effect  is  evidently  dominating  the  effect  of  entrainment  in 
the  Mount  St.  Helens  cloud  because  of  the  large  area  providing  the  source  of 
buoyancy. 

Finally,  we  note  that  the  mosaic  photographs  (Fig.  1)  show  a  prominent 
concentric  ring  structure  in  the  lower  part  of  the  cloud  and  prominent  radial 
protrusions  in  the  uppermost  outwardly  expanding  parts  of  the  cloud.  The 
radial  protrusions  have  also  been  observed  in  the  eruption  cloud  of  Hekla  in 
1947  (Thorarinsson,  1970).  We  are  not  yet  clear  about  the  significance  of 
these  structures. 


III.  THE  DEPOSITS 

* 

The  deposits  formed  from  the  blast  and  the  eruption  cloud  (Clouds  III  and 
IV)  have  been  extensively  studied  (Moore  and  Sissons,  1981;  Ualtt,  1981;  Waltt 
and  Dzurisin,  1981;  Hobblitt  et  al.,  1981;  Walker  and  McBroome,  1983).  This 
work  allows  estimates  to  be  made  of  the  volumes  of  material  in  the  cloud  and 
the  initial  temperature  of  the  pyroclastic  mixture. 

The  surge  and  flow  facies  of  the  blast  deposit  have  an  estimated  volume 
of  0.18  knr*  (Moore  and  Sisson,  1981;  Moore  and  Rice,  1984)  or  0.1  km^ 

(D.R. E.).  The  associated  air-fall  deposit  has  a  volume  of  0.06  km  (0.03  km 
D. R. E.)  within  90  km  of  the  volcano,  as  deduced  from  the  isopach  map  of  unit 
A3  (Fig.  7)  together  with  data  presented  by  Waitt  and  Dzurisin  (1981).  This 
is  probably  an  underestimate,  as  a  significant  proportion  of  the  fine  ejecta 
injected  to  altitudes  of  20  km  or  more  would  be  blown  hundreds  of  kilometers 
downwind.  In  Fig.  7  we  draw  attention  to  the  observation  that  the  isopach  map 
of  A3  covers  an  area  whose  width  in  the  crosswind  direction  is  approximately 
that  of  the  giant  cloud.  The  maximum  thickness  of  6  cm  does  not  coincide  with 
the  volcanic  vent  of  Mount  St.  Helens,  but  occurs  in  an  area  about  12  km  to 
the  north  of  the  volcano.  This  can  be  compared  with  the  Isopach  map  of  the 
plinian  alrfall  deposit,  which  shows  the  point  source  vent  rather  than  the 

broad  zone  of  A3  (Fig.  5b).  We  also  note  that  the  A3  deposit  is  often  rich  in 

accretionary  lapilli,  which  comprise  up  to  50%  of  the  entire  deposit.  In 

addition  to  typical  accretionary  lapilli,  clusters  of  ash  grains  occur  at 

scales  smaller  than  those  normally  considered  ar  accretionary  lapilli.  Liquid 
water,  therefore,  probably  comprised  at  least  10%  by  weight  of  the  deposit. 

The  blast  deposit  consists  of  about  50%  poorly  vesicular  juvenile  dacite 
blocks  and  assorted  lithic  material  of  the  old  volcanic  structure.  Moore  and 
Sissons  (1981)  estimate  initial  temperatures  of  at  least  350°C  for  this 
mixture,  based  on  the  extent  of  charring  of  the  timber.  This  is  consistent 
with  an  equal  mix  of  cold  lithics  and  juvenile  dacite  clasts  (at  850°C)  which 
would  have  a  mean  temperature  of  about  400*0.  Hoblitt  et  al.  (1981)  measured 
temperatures  up  to  277°C  in  the  thickest  pyroclastic  flow  deposits,  but 
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Fig.  7.  (a)  Isopach  aap  of  the  air-fall  deposit  from  the  giant  cloud, 

slapllfied  from  Sisson  (1982).  Contours. are  in  millimeters. 
The  dashed  line  shows  the  outline  of  the  giant  cloud  at  15.55 
GMT.  (b)  Isopach  map  of  proximal  plinian  deposit,  after  Waitt 
(1981),  showing  relationship  to  point-source  vent.  The  outer 
limit  of  the  giant  ash  cloud  is  shown  for  comparison. 


generally  found  temperatures  below  160°C.  Much  of  the  ejecta,  occurring  as 
thin  layers  mantling  topography,  can  be  assumed  to  have  been  emplaced  at  even 
lower  temperatures.  Some  of  the  temperature  decrease  observed  can  be 
attributed  to  the  decompression  of  the  flow  mixture  during  transport  (Kieffer, 
1981).  Conversion  of  initial  heat  and  energy  stored  in  compression  in  the 
gases  into  kinetic  energy  resulted  in  the  high  flow  velocities,  and  some 
cooling  of  the  flow.  Some  of  the  heat  was  also  transferred  to  the  atmosphere 
during  emplacement  as  the  flow  mixed  with  its  surroundings,  and  this  thermal 
energy  source  provides  the  main  cause  of  cloud  ascent. 

Walker  and  McBroome  (1983)  have  drawn  attention  to  the  f ines-depleted 
character  of  the  blast  deposits.  In  particular  the  basal  coarse-grained  layer 
A1  is  strongly  f ines-depleted,  and  they  argue  that  the  ash-fall  deposit  A3 
represents  the  complementary  lost  fines.  The  layer  A2  deposit  is  also  poorer 
in  fines  than  most  flow  deposits.  Irrespective  of  the  relative  merits  of  the 
flow  and  surge  views,  the  grain-size  data  indicate  extensive  mixing  of  air 
with  the  solid  ejecta,  creating  f ines-depleted  flow  (surge?)  deposits  and 
f lnes-enricned  ash-fall  deposits. 


IV.  THEORETICAL  CONSIDERATIONS 


The  observations  in  the  preceding  section  suggest  that  the  formation  of 
the  main  cloud  was  an  abrupt  event  at  the  end  of  blast  emplacement.  This  is 
interpreted  as  a  consequence  of  the  flow  becoming  lighter  than  the  overlying 
air.  Three  main  processes  that  caused  the  flow  to  decrease  in  density  as  it 
moved  away  from  the  source  can  be  identified:  decompression  of  the  blast 
mixture  of  solids  and  gas,  mixing  and  heating  of  the  surrounding  air  into  the 
flow,  and  sedimentation  from  the  flow. 

There  is  general  agreement  that  the  blast  was  triggered  by  the  landslide, 
which  cause  a  major  pressure  drop  within  the  volcano.  Photographic  evidence 
indicates  that  about  500  m  of  rock  overburden  was  removed  (Kieffer,  1981; 
Eichelberger  and  Hayes,  1982).  However,  there  are  two  different  models  for 
the  explosions  that  supplied  the  blast  flow.  In  the  hydrothermal  model  the 
reservoir  is  assumed  to  be  a  mixture  of  water  and  rock  at  or  below  the  satura¬ 
tion  temperature  and  pressure.  In  this  case  the  initial  pressure  may  be 
considered  to  be  between  the  lithostatic  pressure  of  500  m  of  overburden 
(approximately  125  bars)  and  the  hydrostatic  pressure  of  the  hydrothermal 
fluid,  which  might  be  considerably  lower.  In  the  magmatic  model  the  reservoir 
is  assumed  to  be  high-temperature  magma  containing  a  small  fraction  of  ex- 
solved  water  (Eichelberger  and  Hayes,  1982).  The  amount  of  cooling  of  the 
mixture  and  the  change  in  bulk  density  due  to  decompression  are  quite  sensi¬ 
tive  to  the  initial  conditions  and  model  assumed.  Unfortunately  both  models 
give  plausible  explanations  of  the  observed  high  velocities  of  the  blast  flow. 

Hydrothermal  models  require  mass  fractions  of  water  in  the  reservoir  of 
at  least  0.02  to  over  0.1  to  achieve  the  observed  velocities  of  100  to  150 
m/s.  They  also  result  in  substantial  cooling  of  the  blast  flow  mixture  during 
decompression  because  of  the  heat  of  vaporization  and  the  large  mass  fraction 
of  water.  For  example,  in  the  model  of  Kieffer  (1981)  the  mixture  cools  from 
327°C  to  225°C  in  decompressing  from  125  bars  to  atmospheric  pressure. 

Kieffer  argued  that  the  flow  could  become  under  pressured  in  some  regions  and 
therefore  the  cooling  would  be  even  greater.  Her  calculations  suggest  that 


underpresdured  regions  could  approach  the  density  of  the  overlying  atmosphere 
and  that  this  could  cause  the  flow  to  become  buoyant.  In  the  calculations 

O 

presented  by  Kieffer  (1981)  the  flow  density  approaches  2  kg/mJ  in  distal 
regions,  compared  to  an  overlying  atmosphere  with  a  density  of  about  0.8 
kg/ m^. 

In  the  magmatic  model  the  mass  fractions  of  gas  required  to  account  for 
the  observed  velocities  are  much  smaller  and  are  estimated  to  be  less  than 
0.006  by  Elchelberger  and  Hayes  (1982).  As  the  water  vapor  is  initially  at 
high  magmatic  temperatures,  there  is  no  heat  of  vaporization  as  an  energy 
sink,  so  that  the  cooling  on  expansion  is  much  less  than  in  the  hydrothermal 
model.  For  example,  a  mass  fraction  of  water  equal  to  0.006  would  only  cause 
a  temperature  decrease  of  the  bulk  solid/gas  mixture  of  about  13°C  during 
expansion  from  125  bars  to  atmospheric  pressure.  The  main  cause  of  a  decrease 
in  temperature  in  this  model  is  the  incorporation  of  cold  lithic  clasts  during 
the  explosion,  which  could  have  resulted  in  a  temperature  decrease  of  400  to 
450°C  from  the  magmatic  values  (850°C  to  900°C).  A  feature  of  a  magmatic 
model  is  that  the  density  of  the  solid  and  gas  mixture  is  still  high  when 
decompressed  to  one  atmosphere,  because  of  the  small  mass  fraction  of  gas. 

For  example,  at  450°C  and  a  water  content  of  0.006  the  mixture  would  have  a 
bulk  density  of  42  kg/m^. 

Geological  evidence  does  not  favor  a  purely  hydrothermal  model,  as  argued 
in  detail  by  Elchelberger  and  Hayes  (1982).  About  half  the  ejecta  is  juvenile 
vesicular  dacite  and  some  of  the  temperatures  inferred  for  the  flow  appear  too 
high  for  a  shallow  hydrothermal  system.  An  important  magmatic  component  to 
the  explosion  is  indicated.  We  conclude  that  decompression  by  itself  could 
not  have  caused  the  flow  to  have  become  buoyant  and  cannot  account  for  the 
coolings  obsetved  during  flow.  Nevertheless,  a  significant  contribution  of 
hydrothermal  fluids  to  the  explosion  cannot  be  discounted,  and  the  role  of 
decompression  in  reducing  the  flow  density  remains  a  major  uncertainty. 

A  second  process  by  which  the  blast  flow  would  reduce  its  density  would 
be  by  mixing  and  heating  of  surrounding  air  as  the  flow  advanced  over  the 
rough  topography  around  Mount  St.  Helens.  The  complicated  and  pressurized 
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character  of  the  flow  (Kieffer,  1981)  and  the  irregularity  of  the  ground  makes 
impossible  a  theoretical  assessment  of  how  much  entrained  air  would  be 
expected.  As  an  alternative  we  calculate  how  much  air  would  need  to  be 
incorporated  into  the  flow  to  reduce  the  density  to  that  of  the  atmosphere. 

We  make  the  assumption  that  the  mixture  is  homogeneous,  while  recognizing  that 
the  real  flow  is  much  more  complex.  The  idea  of  the  calculations  is  to  assess 
the  average  proportions  of  air  required  to  make  the  flow  buoyant.  The 
conditions  for  this  to  happen  can  be  approximately  estimated  by  considering 
the  heat  balance  between  the  pyroclasts  and  the  air  and  the  definition  of  the 
bulk  density,  which  are  given  respectively  by  the  following  expressions: 


c  M  (9  -0)  =  c  M  (0  -9) 

mme  c  aac  a 


0(1  -  n  ) 
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where  cffl  and  ca  are  the  specific  heat  capacities  of  pyroclasts  (m)  and  air 
(a).  and  are  the  total  masses  of  pyroclasts  and  air  in  the  mixture,  0g 
is  the  pyroclast  temperature,  0  is  the  ambient  air  temperature,  0  is  the 
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mixture  temperature,  a  is  the  ambient  air  density,  0^  is  the  bulk  mixture 
density,  n£  is  the  mass  fraction  of  solids,  and  a  is  the  pyroclast  density. 
With  the  following  definitions, 
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Eqs.  (1)  and  (2)  can  be  considered  to  give 
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Appropriate  choices  of  the  various  parameters  are  listed  in  Table  l.  The 
values  of  air  temperature  and  density  are  estimated  by  using  data  recorded  at 
a  meteorological  station  in  eastern  Washington  on  May  18,  1980.  The  air 
temperature  varied  from  -3°C  to  -23°C  and  the  air  density  from  0.83  to  0.61 
kg/m  over  the  height  interval  2.5  tv>  4  km  above  sea  level.  The  height 
interval  corresponds  to  that  where  mixing  occurred  between  the  blast  flow  and 
the  atmosphere.  The  critical  condition  for  convective  lift-off  is  that  a/8c  = 
1.0.  If  the  total  mass  of  ejecta  is  included  in  the  calculation  (Mm  = 

3.25  x  10^  kg),  the  mass  of  air  at  a/8  =  1.0  is  5.0  x  10^  kg  (625  km^  at  a 

density  of  0.7  kg/nr)  and  the  mean  mixture  temperature  is  152°C.  Over  the 

9 

600-km  area  of  the  blast,  this  mass  of  air  is  approximately  equivalent  to  a 
layer  of  air  1200  m  thick.  These  are  maximum  values  since  it  is  likely  that 
the  flow  became  buoyant  not  only  by  heating  air  but  by  sedimentation.  The 
amount  of  air  required  will  also  decrease  as  the  amount  of  hydrothermal  fluid 
involved  in  the  mixture  increases,  since  the  effect  of  adding  steam  is  to 
reduce  flow  density  further. 

A  minimum  estimate  can  be  obtained  if  only  the  mass  of  ejecta  in  the  air 
fall  deposit  is  used  ( *  0.75  x  10^  kg).  The  calculated  temperature  esti¬ 
mate  remains  the  same,  but  the  air  mass  decreases  to  10^  kg,  equivalent  to  a 
layer  thickness  of  250  m.  This  calculation  is  probably  too  low,  since  we 
suspect  the  volume  estimate  of  A3  is  conservative.  Varying  the  input  parame¬ 
ters  by  amounts  consistent  with  observations  (Table  1)  would  not  change  these 
results  significantly.  The  blast  flow  has  to  mix  with  enough  air  and  sediment 
enough  coarse  ejecta  to  form  a  flow  several  hundred  meters  thick,  to  allow  the 
density  to  fall  sufficiently  to  form  the  major  convective  cloud. 

The  blast  flow  moved  across  the  surface,  progressively  mixing  with  air. 

We  propose  that  the  flow  terminated  when  a  critical  mass  of  air  had  been 
engulfed  and  a  critical  mass  of  sediment  had  been  deposited.  The  blast  flow 
became  buoyant  and  lifted  off  to  form  the  giant  cloud.  The  calculations  are 
broadly  consistent  with  observations.  The  visible  depth  of  the  flow  was  about 
t  km  (Moore  and  Rice,  1984)  and  the  temperatures  of  .he  deposits  were  mostly 
less  than  160°C,  suggesting  that  more  than  half  the  thermal  energy  had  been 
transferred  to  admixed  air. 
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Table  1.  Values  of  Physlc-al  Parameters  Used  In 
Calculations  on  Cloud  Formation 


=  3.25  x  10U  kg 

a  *  0. 8  kg/m^ 

0  =*  A00°C 
e 

0  -  0°C 
a 

a  m  2500  kg/rar1 
cfl  »  993  J  kg-1  K-1 
c_  -  1100  J  kg'*1  K"1 


Thus  we  assume  that  when  the  cloud  became  buoyant  it  had  a  temperature  of 
150°C,  carried  a  mass  of  particles  equivalent  to  that  which  formed  the  air- 
fall  layer  of  7.5  x  10***  kg,  and  had  a  mass  of  air  of  5.0  x  10**  kg.  The 

total  thermal  energy  of  this  material  is  1.07  x  10*  ^  joules,  and  the  density 

<1 

at  an  altitude  of  3  km  (air  density  and  temperature  0.7  kg/nr  and  -10°C, 
respectively)  is  0.48  kg/m'*.  Using  these  estimates  for  the  conditions  of  the 
rising  cloud,  we  can  compare  the  observed  altitude  and  ascent  rate  with 
theoretical  relationships  for  gravitational  convection. 

Morton  et  al.  (1956)  considered  the  instantaneous  release  of  a  source  of 
buoyancy.  For  the  case  of  a  thermal  they  estimate  that  the  final  height  H  is 
a  function  of  the  total  thermal  energy,  Q: 

H  -  KQl/4  (5) 

This  relationship  was  derived  by  considering  the  buoyant  rise  in  a  stably 
stratified  atmosphere,  and  the  constant  K  is  a  parameter  depending  on  the 
environmental  temperature  gradient  and  lapse  rate.  For  a  standard  atmosphere, 
K  «  1.4.  Strictly  this  formulation  should  not  be  valid  for  columns  much 
higher  than  the  stratosphere.  However,  Sparks  (1986)  has  carried  out  numeri¬ 
cal  calculations  on  plume  ascent  which  do  not  strongly  depart  from  predictions 
based  on  Morton  et  al.'s  treatment  of  gravitational  convection,  even  for 
heights  of  30  km.  For  Q  *  1.07  x  10*^  joules,  H  is  25.3  km,  which  is  in  re¬ 
markably  good  agreement  with  the  height  reached  by  Cloud  III  (Fig.  1). 

If  the  cloud  is  regarded  as  a  maintained  plume  of  short  duration  at  25 
km,  column  calculations  (Sparks,  1986)  indicate  a  power  input  of  2.0  x  10*4 
joules/sec  lasting  over  a  9-min  period.  This  is  also  consistent  with 
observations,  since  the  major  period  of  cloud  growth  lasted  10  min. 

The  ascent  profile  for  Cloud  III  can  also  be  interpreted  in  terms  of  a 
simple  physical  model.  Wilson  and  Self  (1980)  analyzed  the  ascent  of  some 
small  eruption  clouds  from  the  volcano  Fuego  in  Guatemala  in  terms  of  the 
ascent  of  a  buoyant  cylinder  of  length  L.  They  equated  the  product  of  mass 
and  acceleration  to  the  total  forces  acting  on  the  cloud  (atmospheric  drag, 
gravity,  and  buoyancy)  to  derive  the  following  expression: 


where  is  the  drag  coefficient,  u  is  the  velocity  of  the  cloud  front,  L  is 
taken  as  the  cloud  height,  g  is  the  acceleration  due  to  gravity,  and  du/dt  is 
the  acceleration  of  the  cloud.  We  have  used  this  relationship  to  estimate  the 

relative  density  g^/a.  The  results  are  listed  in  Table  2  for  Cy  *  1.00  up  to 

18  km  high.  Above  18  km  the  cloud  has  expanded  so  much  that  the  analysis  is 
unlikely  to  be  valid.  The  most  significant  feature  of  the  calculations  is 
that  the  cloud  is  still  buoyant  at  8  km  altitude,  but  has  become  denser  than 

the  surrounding  atmosphere  by  13.3  km.  The  level  of  neutral  density  (8c/a  ■ 

1.0)  is  about  10  km  high.  This  agrees  well  with  radar  data  (Sarna-Woj cicki  et 
al.,  1981)  that  show  that  the  base  of  the  mushroom  cloud  was  at  10  km, 
confirming  that  the  cloud  had  reached  its  neutral-density  level  at  this 
height. 

The  ascent  of  the  cloud  to  25  km  high  is  entirely  due  to  the  column’s 
inertia.  The  cross-section  of  the  column  shows  a  bulge  at  the  top  as  the 
cloud  overshoots  to  25  km  due  to  momentum  and  then  flows  radially  outwards  and 
slightly  downwards  to  form  a  forced  intrusion  between  10  and  20  km  altitude. 
The  anvil  shape  of  the  mushroom  is  characteristic  of  experimental  force  intru¬ 
sions  formed  by  the  ascent  of  buoyant  plumes  into  a  stratified  environment 
(Turner,  1979). 

The  radial  expansion  of  the  cloud  can  be  explained  by  continuity  con¬ 
siderations.  The  mean  radial  velocity,  V^,  can  be  approximately  related  to 
the  total  mass  flux  of  air  and  pyroclasts  at  the  base  of  the  mushroom  cloud 

(at  10  km  high),  where  the  ascending  material  has  the  same  density,  o  ,  as  the 

a 

surrounding  atmosphere,  by  the  following  expression: 

a  u  A 

%  "  -  (7) 

a  2  it  HR 

where  R  is  the  radial  distance,  A  is  the  cross-sectional  area  of  the  ascending 
cloud,  u  is  its  velocity,  a  is  the  mean  air  density  between  the  top  and  base 
of  the  mushroom  cloud  (10  to  20  km  in  this  case),  and  H  is  the  thickness  of 


Table  2 

I 


Values  of  Parameters  Used  in  Calculations  with  Eq.  (6) 


L  (km)* 


g/a 


h  (km) 


u  (m/s) 


du/dt  (m/s^) 


2  2 

the  mushroom  cloud.  A  value  of  A  equal  to  500  km  rather  than  the  600  km  of 

the  blast  area  is  estimated,  since  the  cloud's  cross-sectional  area  decreases 

in  height  up  to  the  height  of  neutral  density  at  10  km.  Taking  values  of  a  * 

0.2  kg/m^,  a  =0.41  kg/nr*,  u=80  m/s,  H  =  10  km,  and  A  =  500  km^,  Eq.  (7) 

<1 

simplifies  to 


1.  3  x  10^ 


R 


(8) 


Figure  8  compares  this  simple  relationship  with  the  data  of  radial  velo¬ 
city  versus  distance,  calculated  by  taking  the  average  velocity  between  each 
5-min  contour  in  Fig.  5.  Although  the  intrusion  of  the  mushroom  cloud  into 
the  atmosphere  is  likely  to  involve  mixing,  the  simple  analysis  appears  to  be 
sufficient  to  model  the  gross  features  of  the  formation  of  the  umbrella- shaped 
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V.  DISCUSSION 


Large  eruption  columns  have  usually  been  associated  with  vertical  explo¬ 
sions  arising  directly  from  volcanic  vents.  The  initial  events  of  May  18, 
however,  have  revealed  a  quite  different  mechanism  involving  the  spread  of  an 
energetic  blast  flow.  The  flow  eventually  mixed  with  enough  air  to  provide  a 
source  of  buoyant  hot  air  and  ash  over  a  600-km^  area.  The  transition  from 
horizontal  density  flow  to  convective  ascent  was  abrupt,  lasting  only  1.5  min. 

There  are  several  features  that  distinguish  eruption  clouds  foiled  in 
this  way.  First,  the  center  of  cloud  ascent  can  be  displaced  away  from  the 
volcano,  by  12  to  13  km  in  this  case  (Figs.  5  and  7).  The  source  of  buoyant 
air  and  ash  is  a  large  area  rather  than  a  point  source,  and  consequently  the 
thickness  of  the  air-fall  deposit  should  be  uniform  in  proximal  areas,  with 
closure  around  the  center  of  the  blast  rather  than  around  the  volcano.  This 
feature  is  observed  with  the  air-fall  deposit  of  the  Mount  St.  Helens  cloud 
(Fig.  7a).  There  is  no  momentum-dominated  basal  region  to  the  column,  as  is 
the  case  for  eruption  columns  generated  by  explosions  (Sparks  and  Wilson, 
1976).  The  cloud  is  only  formed  when  the  flow  material  becomes  buoyant,  and 
must  involve  an  initial  stage  of  acceleration,  rather  than  deceleration  as  in 
explos ion-generated  clouds.  The  Mount  St.  Helens  cloud  only  entrained  fine¬ 
grained  material:  Hoblitt  et  al.  (1981)  and  Walker  and  McBroome  (1983)  report 
only  fine-grained  ash  less  than  2  mm  in  diameter  within  the  deposit.  Coarse 
pyroclasts  were  segregated  into  the  flow  deposit.  The  vertical  velocities  of 
110  m/s  at  8  km  are  sufficient  to  carry  fragments  of  10  cm  diameter  or  more, 
but  the  air-fall  deposit  only  contains  fine-grained  ash  less  than  0.2  in 
diameter.  This  observation  indicates  that  the  coarse  ground  material  is 
segregated  at  low  altitudes,  well  below  8  km  where  the  maximum  velocity  is 
observed. 

Sparks  and  Walker  (1977)  have  proposed  that  pyroclastic  flows  are  associ¬ 
ated  with  the  formation  of  voluminous  air-fall  deposits  as  well  as  ignim- 
brite.  Three  sources  of  fine  ash  for  these  deposits  were  proposed:  (1)  the 
escape  of  fine  ash  and  gas  above  the  collapsing  eruption  column;  (2)  the 
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escape  of  fine  ash  as  a  dense  flow  segregates  near  the  vent;  and  (3)  the 
elutriation  of  fine  ash  from  the  moving  fluidized  flow.  The  behavior  of  the 
May  18  flow  suggests  the  possibility  that  a  substantial  part  of  co-ignimbrite 
ash  clouds  can  be  formed  abruptly  over  the  whole  area  of  the  ignimbrite  during 
the  final  stages  of  emplacement.  The  fine  grain  size  and  wide  distribution  of 
co-ignimbrite  ash  layers  are  more  consistent  with  this  mechanism. 

The  dramatic  expansion  of  the  mushroom  cloud  is  of  general  signifi¬ 
cance.  The  reason  for  the  upwind  flow  and  large  radial  velocities  is  the 
large  total  volume  flux  of  air  and  ash  from  the  vertical  plume  that  supplies 
the  cloud.  The  volume  flux  is  in  turn  related  to  the  temperature  and  mass 
flux  of  solids  (Wilson  et  al.,  1978;  Sparks,  1986).  Irrespective  of  whether 
the  cloud  is  fed  from  a  point  source  (as  in  pllnlan  eruptions)  or  from  a 
large  area  (as  in  the  present  case),  the  development  of  a  large  umbrella  cloud 
should  also  be  a  feature  of  all  eruption  columns  supplied  by  large  mass 
fluxes.  Sparks  (1986)  has  calculated  that  plinian  eruption  columns  with 
heights  greater  than  20  km  and  mass  eruption  rates  greater  than  5  x  10^  m^/s 
should  show  large  mushroom  clouds  and  significant  upwind  flow. 

The  development  of  large  radial  flow  velocities  over  large  horizontal 
distances  has  major  implications  for  the  dispersal  and  grain-size  character¬ 
istics  of  air-fall  deposits.  Carey  and  Sparks  (1986)  have  calculated  the 
trajectories  of  clasts  through  eruption  columns  that  have  large  enough  mass- 
discharge  rates  to  form  a  mushroom  cloud.  The  distribution  of  maximum  frag¬ 
ment  size  is  strongly  influenced  by  transport  within  the  horizontally  expand¬ 
ing  mushroom  cloud.  The  areas  enclosed  by  maximum-grain-size  isopleths  are 
not  influenced  by  either  the  wind  or  the  ascent  in  the  plume.  Maximum-grain- 
size  data  are  controlled  by  the  dimensions  and  velocities  of  the  mushroom 
cloud.  Hence  the  column  height  and  mass  discharge  rate  can  be  estimated  quite 
accurately  from  maximum-grain-size  maps. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "architect -engineer’  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  -xpertise  and  its 
ability  to  stay  current  with  new  deve 1 opmen ' s .  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  wLth 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  die 
research  effort  are  these  individual  laboratories: 

Aerophvsics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  rhemistry,  trace  detection; 
spacecraft  structural  mechanics,  con c ami nat ion ,  thermal  and  strurtuc  >1 
controL;  high  temperature  thermonechan  i.*s ,  gas  kinetics  and  radiition;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  be.im  control,  atmospheric  propag.it  ion ,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  rhemica!  reactions, 
atmospheric  optics,  light  scattering,  state -spec  1 r it:  chemical  reictions  and 
radiative  signatures  or  missile  plumes,  sensor  out  -ot  *! le Id  -ol  -v lew  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  eMect-.  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  transition, 
performance-sensitive  system  design,  distributed  ar.ni lectures  for  spacehorne 
computers,  taul  t -tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  app  L  icat  ions  ,  con.mumc  it  ion  protocols,  and  computer  im.  »tity. 

Electronics  Research  Labor  it  orv :  Microelectronics,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  c  >mmuni cat  ions ; 
microwave  semiconductor  devices,  microwave/mi  1 1 imeter  wave  measurements, 
diagnostics  and  radiomecrv,  microwave/mi  1 l imeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Mater ia  Is  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fricture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
env i ronment  s . 

Space  Sciences  Laboratory:  Magne  t  osphe  r  ic  ,  auroral  and  cosmic  r.;v 
physics,  wave -particle  interactions,  magnetospher ic  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomv, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  e 1 ec t rumagne t l c  and  particulate  radiations  on  space  systems;  space 
inst rumen  tat  ion. 


